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Abstract. We consider a harmonic chain in contact with thermal reser- 
voirs at different temperatures and subject to bulk noises of different 
types: velocity flips or self-consistent reservoirs. While both systems 
have the same covariances in the nonequilibrium stationary state (NESS) 
the measures are very different. We study hydrodynamical scaling, large 
deviations, fluctuations, and long range correlations in both systems. 
Some of our results extend to higher dimensions. 



1. Introduction 

We consider a harmonic chain consisting of N oscillators. The positions 
and momenta of the oscillators are denoted by qj and pj respectively. The 
Hamiltonian is given by 

N+i 2 (a- -a) 2 

3=0 3=0 \i-j\=l, 

i,je{o,...,N+i} 

where we have set the mass of each particle and the nearest neighbor coupling 
equal to 1. We impose boundary conditions qo = qN+i = 0, po = Pn+i = 0. 
The strength of the pinning potential is regulated by the parameter v > 0. 
The energy of site j is given by 

(1-2) ^f + ^f+I £ fe-Si) 2 " 

i: \i-j\=l 

In the absence of any stochastic noise the generator of the dynamics is 
given by the Liouville operator A 

N 

(1.3) A=Y J {pA 3 -l{v 2 -^)<l] 3 dv ] } 

3=1 

where A is the discrete Laplacian with Dirichlet boundary conditions. 

We shall now consider two ways of adding noise to the system. The first 
case we consider is the so-called self-consistent model |11| I12j . Each site 
j £ {!,..., N} is connected to a Langevin reservoir at temperature Tj > 0. 
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The generator of the Langevin reservoir at temperature Tj acting on the 
particle at j is given by 

(1.4) B jjT Tjd 2 p . Pj d Pj . 

The temperatures of the reservoirs at the boundary sites are fixed by the 
conditions T\ = Xg, T/v = T r , while the temperatures of the interior site 
reservoirs are determined self-consistently, by requiring that in the non- 
equilibrium stationary state (NESS) there is no net flux of energy between 
the system and the interior reservoirs. 

The second model we consider is one in which only particle 1 and particle 
N are in contact with heat baths. In the bulk we add a flip dynamics which 
consists of reversing the velocity of each particle at random independent 
Poissonian times. These flips are energy conserving and we call the model 
the velocity flip model. 

In both situations the generator of the total dynamics can be written as 

(1.5) C = A + ~fS + Bi, Ti +B N ,T r 

where 7 > is the intensity of the bulk noise generator S. The phase space 
is f2jv = x Mr and a configuration ((ft, • • • , <?jv>Pi 5 • • • ,Pn) is denoted by 
uj = (q, p). The configuration at time t of the process is denoted by u>(t). 
For the self-consistent chain S has the following expression 

N-l 

3=2 

It can be shown that for any given temperatures T := {Tj > 0; j = 
1, . . . ,N} there exists a unique steady state /x S) g (which depends on the 
TVs) for the dynamics generated by C. Moreover, it can be proved that 
there exists a unique T := T sc such that the self-consistency condition 
lMi,i*~ (Pj) = Tj for j = 2, . . . , N — 1 is satisfied. (If u > 0, the proof is given 
in [11]; the case v = is treated later, in Section l5.2i ) In the sequel we 
fix the temperatures to follow the self-consistent profile T and denote the 
corresponding stationary state by fi^ s or (-) S;S . 

For the velocity flip model the generator S is given by 

1 N 

(S/)(«) = -£[/(^) -/(«)] 
i=i 

where uj 3 is the configuration obtained from oj = (q, p) by flipping pj to —pj . 
Following the lines of the proof given in [3] or [5] one can prove that the 
velocity flip model has a unique non-equilibrium stationary state denoted 

b y Vs,v or Os,v 

The outline of the rest of the paper is as follows. In Section [2] we inves- 
tigate the similarities between the NESS of the self-consistent model and 
the NESS of the velocity-flip model and we show that while the two prob- 
ability measures have the same covariances, (-) SjS is Gaussian but (-) s ,v is a 
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non-trivial mixture of Gaussian states. In Section [3] we obtain the hydro- 
dynamical equations for the velocity flip model. In Section [4] we investigate 
the large deviation fuction and the stationary fluctuations of the energy. 
In particular, we show that long-range correlations are present in the state 
(■}s v This is very different from the NESS of the self-consistent chain: we 
show in Section [5] that long-range correlations are not present in the state 
(■)ss- Some technical matters are deferred to the Appendices. A summary 
of our results is given in [3] . 

2. The steady states 

If the temperatures Tg and T r are equal to a common value T, the steady 
state of the self-consistent chain and the steady state of the velocity flip 
model are both equal to the Gibbs state with temperature T that we denote 
by //jf or (A( ec i> T ). This is a Gaussian measure with covariance C eq (T) = 
TC eq (l). 

In contrast, for Tg ^ T r , the NESS of the two models are different. How- 
ever, it is easy to see that both stationary states are centered, and by the 
results derived in |15| . we know that the two point correlation functions of 
both models coincide. To summarize the result, let us introduce the nota- 
tion (A; B; ■ ■ •) to denote the cumulants of random variables A, B, . . .; for 
instance, (A;B) = {AB) - (A)(B). 

Proposition 1 ([15]). We have 

(A;B) StV = (A;B) S>S , 

for any two random variables A and B linear in {q±, . . . , qN,Pi, ■ ■ ■ ,Pn}- 

It follows from Proposition [1] that both models satisfy Fourier's law with 
the same value of the conductivity. Indeed, since the energy current 

fx = ~2^ Px +Px-l)(9x ~ Qx-i), x = 2,... ,N- 1, 

is a quadratic function of the momenta and positions, its averages over (-} s s 
and (-) s ,v are equal. This implies that the conductivity of the self-consistent 
model, k S;S (T), and the conductivity of the velocity-flip model, «; SjV (T), both 
defined by 



k S; .(T) = lim lim 



T L ,T r ^-T N^-oo Tg - T r 

are equal to the same value k(T). The existence of these limits have been 
proved in [11] for the pinned (y > 0) self-consistent chain and in [5] for the 
unpinned {u = 0) velocity-flip model. By using Proposition [H we get the 
existence and equality of these limits for both models in the pinned and 
unpinned cases. The value of k(T) is given by 

1/7 



(2.1) «(T) 



2 + V 2 + 0/2(1/2+4)' 
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The steady state (-) S s of the self-consistent chain is Gaussian. The prob- 
ability measure (-)s,v is not Gaussian but we will show it is a mixture of 
Gaussian states. 

Proposition 2. There exists a probability measure p StC ov whose support is 
included in the set £ composed of positive definite symmetric matrices, such 
that {')s iV is given by 



(■)s,v = J G C {-) Ps,cov(dC). 

Proof. See Appendix [A] □ 

There is a simple non-trivial consequence of this property: 

Corollary 1. Let /(q, p) be any function of the form 

N 

/(q>P) = ^{aiQi + biPi) 

i=l 

where (a±, . . . , a^, b±, . . . , &at) are arbitrary real numbers. Then, the follow- 
ing inequality holds: 

(f A )s, v >3[(f 2 ) StV ] 2 = {f 4 ) s , s . 

Proof. If Z € M. d , d > 1, is a centered Gaussian random variable then 
E(.Z 4 ) = 3 [E(Z 2 )] 2 . Thus, for any Gaussian measure Gc with covariance 

matrix C, we have Gc{f A ) = 3 [Gc(f 2 )] 2 ■ By Proposition [2] and Jensen's 
inequality, the result follows. □ 

3. Hydrodynamical scaling limit of the velocity flip model 

We have to distinguish two cases according to whether v = (unpinned) 
or v > (pinned). The unpinned case is similar to that investigated by one 
of the authors in [2] . 

3.1. The unpinned chain. When v = it is useful to define the more 
convenient coordinates r x , < x < N by r x = q x +i — q x , and set rjy + i = 
0. The phase space is then identified (with some abuse of notation) as 
Otv = M^" 1 " 1 x 1S> N and a configuration (ro, . . . , r^,p±, . . . ,pn) is denoted by 
ijj = (r, p). The configuration at time t of the process is denoted by ui(t). 
In terms of these new variables the Liouville operator can be written as 

N N 
x=0 x=l 

The bulk dynamics conserves two quantities. The first one is the energy 
%. The second one is the deformation ^2 x r x of the lattice. Because of 
the Langevin baths the total energy is no longer conserved by the total 
dynamics; the situation is different for the total deformation ^2 X=0 r x since 
our boundary conditions go = Qn+i = fix it to at any time. Nevertheless, 
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the field r(t) = {r x (t) ; < x < N} fluctuates, and this second conservation 
law has to be taken into account in the hydrodynamic analysis. Observe 
that the conservation of total momentum is no longer valid because it is 
broken by the noise. 

The energy at a site x G {1, . . . , N} is now given by 

(3.1) ^ = 4+4+ r ' _1 



2 4 4 



and £q = jTq , £n+i = j^n ■ The local conservation of the two conserved 
quantities is expressed by the equations 

£ x (t)-£ x (0) = f ds {f x (s) - f x+1 (s)} ds 
(3.2) Jo t 

r x (t)-r x (0)= / ds {f x (s)-f x+1 (s)}ds 



where the current of the energy j x and of the current of the deformation of 
the lattice j x are given by 

(3-3) H = --r x _i{p x +p x -i), j x = -Vx ■ 

The equilibrium Gibbs measures of the infinite system are parametrized 
by the inverse temperature (3 = T~ l and the pressure (or tension) r. They 
are given by 

JV+l 

(3.4) fi TtT (oj) = Z(T, r)- 1 11 exp{-P(£ x - rr x )} 

2=0 

and sometimes denoted by {-)t,t- We note that the following relations hold 
for x e {1,...,JV} 

(3.5) (pl) TjT = T, (£ x ) t , t = T + t 2 /2, {r x ) T>T = T. 

There is a second natural parameterization of the equilibrium measures 
by the averages of the two conserved quantities. Let E, R be the functions 
defined by 

(3.6) E(T, t) = T + r 2 /2, R(T,t) = t. 

For any R G M, E > ^R 2 , we denote by fi^ r the canonical equilibrium 
measure ij,t,t with r, T such that E(T, r) = E, R(T, r) = R. Explicitly, 
these are given by r = R and T = E — ^R 2 > 0. 

To establish the hydrodynamic limits corresponding to the two conserva- 
tion laws, we look at the process with generator N 2 C, i.e., in the diffusive 
scale. Assume that initially the process is started with a Gibbs local equilib- 
rium measure Aeo (•)>«<)(•) associated with a macroscopic deformation profile 
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Uo(q) and a macroscopic energy profile £o(<?) : 
( 3 - 7 ) Ae (-),«o(-) 



1 N+l 



[ exp{-/3 (^/^)(^x - T (x/N)r x )}dr x dp x 



where To = /Sq -1 and ro are the temperature and tension profiles correspond- 
ing to the given energy and deformation profiles through the correspondence 
defined in (|3.6p . The corresponding expectation is denoted by (•). Both pro- 
files are assumed to be continuous and then we have for any macroscopic 
point q £ [0, 1] 

(3.8) Jim (r [Ar j(0)) =u (q), lim (£ [Nq] (0)) = e (q) . 

We show that at any later (macroscopic) time t 

(3.9) lim (r [Nq] (N 2 t)} = u(q,t), lim (S [Nq] (N 2 t)) = e(q,t) , 

TV— >oo 1 ' TV— >oo 1 J 

where u, e are solutions of the following macroscopic diffusion equation 
(3.10) 



dtu = 7 1 d q u 



d t e={2 1 y 1 dl(e + u>/2) 

with the initial conditions n(-,0) = uq(-), e(-,0) = so(-) and boundary con- 
ditions 

(d qU ) (o,f) = (a q u)(i,t) = o, 

(3.11) / u 2 \ . / u 2 



{0,t)=T t , U- — )(l,t)=T r . 

Taking the limit t — > oo in these equations we obtain the typical macro- 
scopic profiles of the system in the nonequilibrium stationary state i.e., 
a flat deformation profile u = and a linear profile T interpolating between 
T £ and T r , 

(3.12) f(q) = T e + (T r - T t )q , 

for the energy profile. 

3.1.1. Derivation of Macroscopic Equations. We sketch a proof of (|3.10p at 
a somewhat informal level. It is based on the following two part analysis. 
A fully rigorous treatment would involve dealing with several non-trivial 
technical problems. We will not do this here but see [TJ [3] for a rigorous 
proof in a similar context. 

(1) Local thermal equilibrium (LTE): Local thermal equilibrium corre- 
sponds to each given small macroscopic region of the system being 
in equilibrium, but different regions may be in different equilibrium 
states, corresponding to different values of the parameters (see (13.71) ). 
Let us consider a very large number, say 2L+1, of sites in microscopic 
units (L S> 1), but still an infinitesimal number at the macroscopic 
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level ((2L + 1)/N < 1). We choose hence L = SN where S < 1. 
We consider the system in the box A^(x) composed of the sites 
labeled by x — L, . . . ,x + L. The time evolution of the 2L + 1 oscilla- 
tors is essentially given by the bulk dynamics since the variations of 
deformation and energy in the volume containing the 2L + 1 oscil- 
lators changes only via boundary fluxes. After N 2 microscopic time 
units, the system composed of the 2L + 1 oscillators has relaxed to 
the micro-canonical state ^e q (tN 2 ),f q (tN 2 ) at q = x/N corresponding 
to the local empirical deformation r q (tN 2 ) and the local empirical 
energy e q {tN 2 ) in the box Al(x). Indeed, we can divide the observ- 
ables into two classes, according to their relaxation times: the fast 
observables, which relax to equilibrium values on a time scale much 
shorter than N 2 and will not have any effect on the hydrodynamical 
scales, and the slow observables which are locally conserved by the 
dynamics and need much longer times to relax. Let (j){u)) be a local 
observable and let <f> x = r x <j> denote the same observable translated 
by x. Thermal local equilibrium corresponds to {4>^ q ](N 2 t)} being 
close to its value in ^e q (tN 2 ),f q (tN 2 )- I n the thermodynamic limit we 
take first N — > oo and then 5 — > (we recall that 5 is related to 
the size of the box A^(x) by L = 5N), we have, by equivalence of 
ensembles, 

(<t>{Nq]{N 2 t)) W h(q,t),u{q,t){4>) 

since 

e q {tN 2 ) W£(},i), f q (tN 2 ) f*u(q,t). 

(2) First order corrections to LTE: Since we observe the system on a 
diffusive scale, local thermal equilibrium alone is not sufficient to 
determine the form of the hydrodynamic equations. The second 
key ingredient is the computation of first order corrections to LTE. 
These corrections are expressed (in the bulk) by what is called a 
"fluctuation-dissipation equation" in the hydrodynamic limit math- 
ematical literature. In general, such equations are not explicit, but 
this happens for our system. We begin by noting that the currents 
are given, on the microscopic level, by 

3% = -{V<l>) x + C{h x ), f x = - 1 -\Vr) x - 1 +£{ 1 - l p x ), 

) 10 1 

4>x = (27)" {Px-i + r x -ir x -2), K = -7" 3% , 
where C is defined in (II. 5|) and (V(f>) x = (p x +i — (f> x . 
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To get the hydrodynamical equations, we let G : [0, 1] — > K be a test 
function. Then we have 

N N 

N- 1 £ G(x/N)(£ x (tN 2 )) - N- 1 £ G(x/N)(£ x (0)) 

x=l x=l 
N ft N 2 



^ 2 E / G'(x/N)(f x (s))ds 

N ft N 2 

N- 3 Y] / G"(x/N)(Ms))ds 

x=l J ° 

N ft N 2 

+N- 2 J2 G'(x/N)({Ch x )(s))ds 

x=l J ° 

N- 1 V / G"(x/N)(MsN 2 ))ds 

x=l J 

N 

+N~ 2 Y J G'{x/N) {(h x (tN 2 )) - <M0)>} 



x=l 
N r-t 

w iV _1 V / G // (x/iV)((/>x(s^ 2 ))(is + 0(^ 1 ). 

To get the first equality, we have used local conservation of energy (|3,2p and 
a discrete integration by parts. To get the second one, we used (|3.13p and 
a second discrete integration by parts. To get the third equality, we used 
Lemma [2] in Appendix B. 

By LTE the term ((p x (sN 2 )) is asymptotically equivalent in the thermo- 
dynamic limit to 

P>e(q,s),r{q,s)(<i>x) 

where q = x/N. A simple computation shows that 

MT.rfa*) = (27)" 1 (r 2 + T) = (27)- 1 [E(T, r) + R 2 (T, r)/2] 

so that 

AiUifa*) = (2-y)- 1 {E + R 2 /2) . 

This closes the evolution equation for the energy field. 

A similar computation can be performed for the deformation field giving 
the form of the hydrodynamic equations A3. 10f> in the bulk. These equations 
are, as already noted, to be solved subject to the boundary conditions (I3.11|) . 
The reason for this is that the time relaxation of the Langevin baths is of 
order one and we look at the process at time tN 2 . On this time scale the 
system has reached LTE with temperature Ti at the left and temperature 
T r at the right. Since the average of the kinetic energy at (macroscopic) 
time t and position q G [0, 1] is given by (e — u 2 /2){q,t), we get the second 
equation in (13. lip . The Neumann boundary conditions for the u-field are a 
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consequence of the conservation of the total length of the chain ^2 X=0 r x (t) = 
qN + i(t) — qo(t) which remains equal to during the time evolution. Thus 
we have Jj ^ u(q, t)dq = for any time t. Since in the bulk u is a solution of 
the heat equation, by taking the time derivative of the previous equality we 
get f(t) := d q u(0,t) = d q u(l,t) (the fluxes at the boundaries compensate). 
We have to show that f(t) = 0. Observe that 

(3.14) £pi = (n - r ) - (7 + l)pi, Cr = p 1 . 
Thus, by Lemma [2] in Appendix [Bj we have for any t > 0, 

^{( Pl (tN 2 )) - ( Pl (0))} = jf* N [(n( S N 2 )) - (r (sN 2 ))] ds 

-( 7 + l)JV / ( Pl (sN 2 ))ds, 
Jo 

j\ Pl (sN 2 ))ds = -L{( ro (tN 2 )) - <r (0)>} 

The second equality shows that N fQ(pi(sN 2 )}ds vanishes as N goes to 
infinity. The left hand side of the first equality is of order 1/N and since 
fo N [ri(sN 2 ) — ro(sN 2 )] ds converges to f Q (d q u)(0, s)ds = f Q f(s)ds, we 
have f(s) = for any s. 

3.2. The pinned chain. When v > 0, the energy current j* satisfying 
M x = 3% - f x+l is given by 

fx = ~2^ Px +Px-i){qx - qx-x) ■ 

Because of the presence of the pinning, the bulk dynamics conserves only one 
quantity: the energy %. It follows that the Gibbs equilibrium measures of 
the infinite system /x^ are fully characterized by the temperature T = 
Under /i^ the p x are independent of the q x and are independent identically 
distributed Gaussian variables of variance T. The q x are distributed accord- 
ing to a centered Gaussian process with covariance UriQxQy) = F(x — y) such 
that 

(3.15) [, 2 -A]r( Z ) = Tl {0} (z). 

In particular, we have fixiPx) = ^t{£x) = T. 

The "fluctuation-dissipation" equation (|3.13p . which gives the first order 
corrections to local thermal equilibrium, is now 

(3.16) 3% = -{V<t>)x + C(h x ) 

with h x = — , y~ 1 j x , as before, but we need to choose now 



(3.17) <j) x = — [pl_! -(v+ - q x q x -2 + Qx-iQx + qx-lQx-2] 



10 C. BERNARDIN, V. KANNAN, J. L. LEBOWITZ, AND J. LUKKARINEN 



Assume that the system is initially distributed according to a Gibbs local 
equilibrium measure associated to the energy profile eo(q), q £ [0,1], and 
define s(q,t) as the evolved profile in the diffusive scale, i.e., 

(3.18) e(q,t) = lim (£ [Nq] (tN 2 )) . 

Then the arguments of Subsection 13. II and (|3. 16[) show that e is the solution 
of the following heat equation 



(3.19) 



d t e = d q (D(e)d q e) , 
e(;0)=e (-), 
{e(0,t)=T e , e(l,t)=T r . 



The diffusivity (in general called the conductivity in the context of heat 
conduction) D(T) is given by 

(3.20) g(r) = - ^y 



By (HTTP and ([37T5]) we get that 

M<t>*) = ^[r-(^ + i)r(o)-r(2) + 2r(i)] 
= A[r -z, 2 r(o)-(Ar)(i)] 
= i[r -^ 2 (r(o) + r(i))] 

and {u 2 + 2)r(0) - 2r(l) = T, T(0) = T J^{u 2 + asm 2 ^))" 1 ^. Thus we 
have 

1/7 

(3.21) D := D(T) 



2 + v 2 + vV(^ 2 + 4) 

Observe that, as expected, the value of D is equal to the value of the con- 
ductivity k(T) in the NESS (see (I21]> ). 

When t goes to infinity e(q, t) converges to the linear profile T(q) = Tg + 
(T r — Tg)q. The latter is the typical profile observed in the stationary state 
("}s v We note finally that, since the self-consistent model does not conserve 
energy in the bulk, we do not expect any autonomous macroscopic equations 
in that model. 



4. Energy Fluctuations in the stationary state of the velocity 

flip model 

4.1. The unpinned chain. In order to obtain the form of the energy fluc- 
tuations in the steady state we first derive the form of the energy dynamical 
fluctuations in the diffusive scale. We start again with the unpinned case. 
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4.1.1. Dynamical fluctuations for the unpinned case. The deformation and 
energy fluctuation fields are defined by 

1 N 

Uf{F) = —J2F(x/N) [r x (tN 2 ) - u(x/N,t)] , 
1 N 

y t N (G) = -=^G(x/iV) [£ x (tN 2 ) - e(x/N,t)] , 



x=l 

where F, G : [0, 1] — > M are smooth test functions vanishing at the bound- 
aries and u,e are the solutions of the hydrodynamic equations (13.10j) . For 
notational convenience we also introduce the temperature profile T(q, t) as- 
sociated to e, u and defined by 

(4.1) r(ff,t) = e (g,f)-*M 

We assume that the process is distributed according to a Gibbs local 
equilibrium state with deformation profile uq and energy profile £q. One 
shows easily that as N goes to infinity (TZq^q) converges to some limit 
field (TZq, 3^o ) but since we do not need the explicit form of the latter in the 
sequel we do not compute it here (see [19]). 

Let a, b, c be the space-time dependent functions given by 



(4.2) 




T 2T 
v 7 ^' V 7 ' 

We claim that as N goes to infinity {K N ,y N ) converges to (TZ, y) given by 
the solution of the following coupled equations 

(4.3) 



d t n = T x dln- d q (cWi; 



q 

d t y = (27)- 1 (3 g 2 (uK) + d 2 q y) - d q (aWi + bw 2 ) , 

where Wi,W 2 are two independent space-time white noises. The initial 
conditions, independent of W±, W 2 , are given by 1Zo and 3^o- 

The proof of (|4.3p is the following. For any Gibbs equilibrium measure 
Ae r with mean energy E and mean deformation R, and for any local ob- 
servable / depending on the configuration uj, let f(E,R) = As^(/)- Using 
the fluctuation dissipation equations (|3.13p , standard stochastic calculus and 
the hydrodynamic equations fj3. 10I) . it is not difficult to show that 
(4.4) 

K?(F) - K»(F) « 7" 1 f* K»{F")ds + M?(t), 
y t N (G)-y?(G) » £ -^J2 x G"(x/N) U x (sN 2 ) - fae{jj ,s), 



N 

+ M»(t) 
where 

(4.5) 4> x (E,R) = (2 1 )- 1 (E + R 2 /2) 



ds 
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and Mi , A4^ are two martingales. By using (IB. 5ft of Appendix |B] and the 
LTE assumption, we get that, as N goes to infinity, 

([M? (tf) -> ^j\sj\q[H'(q)] 2 T(q,s), 

( [M%{t)] 2 )^-f ds C dq[G'(q)] 2 T(q, s) (T(q, s) + 2u\q, s)) , 
7 Jo Jo 

{M N (t)M?(t)) -+ - f ds [ dq[H'(q)G'(q)]T(q,s)u(q,s). 



7. 

This implies that the martingales (Mi , M 2 ) converge as N goes to infin- 
ity to (cWi(H'), (aWi + bW 2 )(G')) with Wi,W 2 two standard independent 
space-time white noise and a,b,c are the functions introduced in (|4.2p . 

The first equation of (|4.4p is closed and the martingale convergence result 
gives the first equation of f(4.3j> . The second one is not and the closure is 
obtained through a nonequilibrium Boltzmann-Gibbs principle |20j which 
generalizes the standard equilibrium Boltzmann-Gibbs principle [13]. Let 
us explain what it means. Observables are divided into two classes: non- 
hydrodynamical and hydrodynamical. The first ones are the non-conserved 
quantities and they fluctuate on a much faster scale than the conserved 
ones. Hence, they should average out and only their projection on the 
hydrodynamical variables should persist in the scaling limit. Thus there 
exist macroscopic space-time dependent functions C(q, s), D(q, s) such that 
for every test function J(q), q £ [0, 1], and any local observable / depending 
on the configuration u), 

1 ft N 

lim — dsJ2 J{x/N) (r x f)(sN 2 ) - f(e(x/N, s), u(x/N, s)) 

-C(x/N, s) (r x - u{x/N, s)) - D(x/N, s) (S x - e(x/N, s))} 

= 

in mean square norm. The functions C and D depend on the macroscopic 
point q = x/N and on the macroscopic time s. In order to compute them, 
we assume local thermal equilibrium. Around the macroscopic point q, the 
system is considered at equilibrium with a fixed value of the deformation 
u(s,q) and of the temperature T(s,q). The functions C,D are then com- 
puted by projecting the function r x f on the deformation and energy fields. 
The values of C and D are given by 

C(q, s) = (d R f) (u(q, s),e(q, s)), D(q, s) = (%/) (u(q, s),e(q, s)). 

Using the Boltzmann-Gibbs principle we see that we can close the time 
evolution equation for the energy fluctuations field. By (14. 5p . we have 

C(q,8) = -^-, D(q,s) = — 
and we obtain the second equation in (|4.3|) . 



HARMONIC CHAINS WITH BULK NOISES 



13 



4.1.2. Stationary fluctuations. We are now ready to compute the energy and 
deformation fluctuations in the steady state, i.e., in the large iV limit of 

1 N 

n»(F) = —Y,F(x/N)r x , 

* x=l 

1 N 

* x=l 

where the configuration (r, p) is distributed according to the steady state 
(■}s,v anci T is t ne temperature profile (|3.12p . 

We show in Appendix ICl that {y ss ,1Z ss ) = Huin-^oo {y^i^ss) are two 
independent Gaussian fields with covariance given by 
(4.6) 

{[T^ss(F)] 2 ) s ,v = /J T(q)F 2 (q)dq 

{[y ss {G)} 2 )^ = £ f\q)G 2 (q)dq + (T r - T t f £ C^XC-Ao)" 1 ^)^ 

where Ao denotes the Laplacian with Dirichlet boundary conditions on [0, 1]. 

4.2. The pinned chain. Our goal is to estimate the probability that in the 
stationary state the empirical energy profile defined by 

N 

(4.7) o N ( q ) = Y,£xi[x/N,(x+i)/N)(q), ?e[o,i] 

2=1 

is close to a prescribed macroscopic energy profile e(q) different from T(q), 
i.e. we want to find the large deviation function (LDF) for the NESS. 

At equilibrium, Ti = T r = T = j3 , the stationary state (-) ss coincides 
with the usual Gibbs equilibrium measure //^ ' cq defined by 

No 1 ( N \ N 

d^' eq = - exp I -p^£ x J Y[ dq x dp x 

' \ x=l / x=l 

where Zjy{T) is the partition function. By the usual large deviations theory 
(see e.g. |14j ) we have that for any given macroscopic energy profile e(-) 

(4.8) ^' cq (9 N (-) ~ £(•)) ~ e~ NV ^ 
where 

e(?) 1 . / e(g) 



(4.9) V cq (e) = J |^-l-lo g 1^1 1 . 

The quantity V cq is the large deviation function (LDF) in the canonical 
ensemble. The LDF coincides for equilibrium systems with the difference 
between the free energy of the system in LTE and the equilibrium free energy. 
It will thus vanish when e(q) = T . 

The purpose of the macroscopic fluctuation theory (MFT) of Bertini et 
al. [6] is to obtain the probability of a large deviation for boundary driven 
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diffusive systems. It predicts that out of equilibrium (Tg ^ T r ) a large 
deviation principle in the following form holds 

(4.10) fi v (*"(•) ~ e(-)) ~ e- NV ^) . 

Moreover, the rate function V can be identified with the so-called quasi- 
potential, a quantity introduced in the context of stochastically perturbed 
dynamical systems [17], and generalized by Bertini et al. to the infinite 
dimensional context; V is then the LDF for the energy profile. 

The explicit form of the quasi-potential is in general unknown. But the 
MFT has been rigorously proven in the context of the boundary driven sim- 
ple exclusion process \W\ [TB] and formally derived for few other systems 
[21 El [9]. The main interest of the MFT is that it does not rely on the 
microscopic properties of the system studied but only on two macroscopic 
quantities: the diffusivity and the mobility [8]. In particular, the LDF V 
depends only on these two data. 

For the system we are interested in, the diffusivity D(T) is given by (|3.2ip . 
By the Einstein relation the mobility x(T) is equal to x(T") = D(T)a(T) 
where cr(T) is the static compressibility defined by 

(4.11) a{T):=Y J ^m-T){8 x -T)) . 

iceZ 

A simple computation shows that cr(T) = T 2 and Theorem 6.5 of [8] applies. 
It follows that V(-) is given by 

(4.12) V ( e) = ^ * - 1 - log (^—j - log 



where F is the unique increasing solution of 

e 



(4.13) 



d 2 n F F 



(d q F) 2 F 2 ' 

F(0) = T t , F(l) = T r 



Thus the function V is independent of the pinning value v 2 and of the 
intensity of the noise 7. In fact, it coincides with the LDF of the KMP 
model considered in 

It is now easy to derive the Gaussian fluctuations of the empirical energy. 
We consider a small perturbation, e = T + 6h, of the stationary profile e. 
The functional V has a minimum at T so that 

V{e) = V(T) + ^S 2 (h, C~ l h) + o{5 2 ) 

The operator C is the covariance for the Gaussian fluctuations of the empir- 
ical energy under the invariant measure ^^ v . The computations are exactly 
the same as in [9], section 5.2, and we get 

(4.14) C = f 2 l + (T r -T e ) 2 (-A )- 1 
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where Ao denotes the Laplacian with Dirichlet boundary conditions on [0, 1]. 
Therefore, it coincides with the expression (14. 6ft even though (|4.6j) is derived 
assuming v = 0. 

4.3. Higher dimensions. It should be possible to extend the previous re- 
sults to higher dimensions, with or without pinning, apart from the deriva- 
tion of the LDF performed in subsection 14.21 A possible difficulty for v = 
is that the Gibbs equilibrium measure has long range correlations in the q 
variables and they are even unbounded in the volume for d = 1,2. At equi- 
librium, these problems can be avoided by considering the corresponding 
gradient fields [18]. For instance, in d = 1 these are given by the variables 
fx = Qx+i — Qx, which have a product equilibrium measure. 

5. Energy correlations in the steady state for the 
self-consistent chain 

In this section we consider the self-consistent chain and mainly follow the 
notations and results in |llj . We show that the energy correlations in the 
steady state are not long range in the sense that 

1 1 N+1 

(5.1) y = jj £ {£ 3 r,£^ T r)+R N , 

j',j=o 

where {Tj c } denotes the self-consistent profile and -Rat is a remainder of 

order 1/yfN in the pinned case and A -1 / 4 In 2 A in the unpinned case. The 
fluctuations are thus dominated by the local equilibrium term, whose pro- 
file is well approximated by a linear profile connecting the fixed boundary 
temperatures 

5.1. The pinned self-consistent chain. We begin with the pinned case 
v > 0, since the necessary estimates in this case have been proven in [llj . 
The local energy at site j G {0, . . . , N + 1} can be written as 

(5.2) £ j = ±pj + ±q r A®q, with 

7V+1 

45 = = X = y] + l[ X = y]Y J B® - B$ , 

y'=o 

B$j = t[\x -y\ = l]i(l[x = j] + l[y = j}) , 

where t[P] = 1, if the condition P is true, and t[P] =0 otherwise. For 
any choice of temperatures of the heat baths T = {Tj ; j = 1, . . . , N} the 
NESS jU s ,s , denoted to simplify notations by (•), is Gaussian with mean 
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zero. Hence we have 

(5.3) ±(£j'£j) 

= (PfP 2 ,) + (q T A^qq T ,4^q) + (pjq^q) + (pfq T A^q) 
= 4(£ r )(£ J )+2(p fPj ) 2 + 2 A^,A^( qx ,q x )(q y ,q y ) 

x'y'xy 



■f,y 



■<7/ 



where the sums go only over {1,2,..., iV}, as = qo = qN+i = Po = PN+i- 
In particular, at thermal equilibrium with all the temperatures Tj equal to 
T, we have 



(5.4) 



2{£ f ;£ j )(^ T ) = T 2 (l[f = j g {1,2, ... ,N}] 



+ Tr 



A^'\-A + u 2 )- 1 A^(-A + 



.2\-l 



Since v > 0, (— A + v 2 ) X y is exponentially decreasing in \x — y\ with an 

exponent independent of N. Therefore, using the fact that A X y is zero 
unless \x — j\, \y — j\ < 1, 



N+l 



(5.5) 



p A«')(-A + v 2 y l A^{-t\ + v 2 )- 1 

j'=0 



is uniformly bounded in N and j. 

Let T = {Tj C ; j = 1, . . . , N} denote the self-consistent temperature pro- 
file, which has been proven to lie between the boundary temperatures Ti 
and T r , and thus is uniformly bounded in N. For convenience, define also 
To = Ti and Tjv+i = T r . Then, 

N+l 

(5.6) £ \(£ f ; £j )^\=0(N). 
j',j=o 

Thus the fluctuations of the energy at the NESS satisfy 

.. 1 N+l N+l 

(5.7) -{%,%) = - J2 = nH (£f,£ 3 ) M) + RN, 

j'J=0 j',j=0 

where the first term is O(l) and, by (|5.3p . the remainder is bounded by 
1 N+l 

(5.8) \r n \ < — £ mm 2 ) + E 4>S¥<^><«>) 



j'd=o 



x'y'xy 
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Here we have used the shorthand notation Sj to denote the difference be- 
tween a quantity evaluated with NESS and the equilibrium distribution at 
temperature Tj, i.e., 

(5.9) Sj(X) = (X) - (X)( cq ' T i\ and 

(5.10) 8j((X)(Y)) = (X)(Y) - (x)( cq ' T ^(y)( cq ' T ^ . 
Since 

(5.11) Sj((X)(Y)) = {Y)5 3 {X) + {X)^6j{Y) , 
we have 

1 N 

(5.12) \R N \< Y. T M^)) 



JV+l 



j',j=0 x'y'xy 
1 N+1 

+ ^a7 E \(PfPj) 6 j((Pj'Pj))+ E A xi A{ ii^'^) 5 ji(%'<ly)) 



j',j=0 x'y'xy 

+ E A i?(PjQy)^((PjQcc)) + Yl A( Jy{Pj"ly} S j((Pf<Ix}) ■ 

xy xy 

It is proven in [11] that Tj is linear up to corrections which are uni- 
formly bounded by TV" 1 / 2 , and that \8j((XY))\ < CiV" 1 / 2 for any choice 
of X, Y G {qj,Pj}j where the constant C depends only on the fixed bound- 
ary temperatures and on v, 7. In particular, it was proven that en '■= 
maxi</ c< jv l^fc+i ~~ Tk\ = 0{N~ 1 / 2 ). In fact, as we shall show next, the pre- 
vious estimates can be straightforwardly strengthened to an upper bound 
with exponential decay in the difference \f — j\. 

We use the results proven in Sec. 2-4 of |llj . There it is shown that 

to each of the pair (X,Y) in (q,q), (q,p), (p,q), (p,p) corresponds a B x k y 
which is independent of the temperatures of the heat baths and using which 
(X x Y y ) = Ef=i B x k y ] T k in the NESS of any temperature profile {T k > 0}. 
In addition, it was shown that 

(5.13) Bg = f N (x -k,y-k) + f N (x + k,y + k) 

- f N (x - k,y + k) - f N (x + k,y - k), 

where the function /jv is exponentially decaying in both arguments in the 
precise sense that there are ^-independent constants C, a > such that 

(5.14) \f N {x,y)\<Ce- a ^' + \y\'\ 

with the periodic distance \x\' := \x mod 2(iV+l)| where x mod 2{N+1) G 
{-N, -N + 1, . . . , N + 1}. For x, r G {1, 2, . . . , N}, one obviously has \x - 
r\ < N — 1 implying \x — r\' = \x — r\. But \x + r\' = x + r only for 
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x + r < N + 1, and else \x + r\' = 2N + 2 — x — r. For x + r > N + 1 thus 
|x + r|' = 2(iV + 1 — x) + x — r = 2(N + 1 — r) + r — x which implies that 
|x + r|' > \x — r\, and this inequality is obviously true for x + r < N + 1. 
Thus we can conclude that for any choice of the signs 

(5.15) 

\x ± r\' + \y ± r\' > \x — r\ + \y — r\ > -(\x — y\ + \x — r| + \y — r|) , 

where the second inequality follows from the triangle inequality. Therefore, 
we have now proven that 



(5.16) 



B^\ < 4Ce~ k a (\ x -y\ + \ x - r \ + \y- r \) . 



'xy 



Using these estimates, we can thus improve Corollary 4.2 of [11] to con- 
clude that for any temperature profile and for any x, y,j in {1, 2, . . . , iV}, 



(5.17) 



N 



(X x Y y ) - {X x Y y )^\ < \B<8\\T k - Tj 

k=l 

N 

< Ac^^ lia{lx ~ vl+lx ~ kl+ly ~ kl ¥ ~ j\^N 

k=l 

< ACeNe-^-y^e-^Hnl + \x- j\) 



net 



where en 
also for j 



- maxi<fc<jv \Tk+i — Tk\. Clearly, the final bound is then valid 
0, N + 1 simply since To = T\ and Tjv+i = T/v. Thus for the 



self-consistent profile there is C > independent of N such that 



Sj((X x Y y ))\ < -^=(1 + \x - j|)e-H^I 



(5.18) 



In addition, (|5.16p implies that there is C" > such that for any T > 0, 



(5.19) 



x-y\ 



k=l 



and for the selfconsistent NESS 

(5.20) \{X x Y y )\ < C" maxiT^e- 1 ^-^ . 

Applying these bounds in (|5.12p shows that there is C independent of N 
such that 
(5.21) 



\Rn\ 



+ 



C r 1 x 



e 2 



l a (\ X '- X \ + \y'-y\) 



N 



j \x-y\<l 



\x-y\<l,\x'-y'\<l 



N 



e 2 



\a(\j'-x\+\j>-y\) 



j' \x-y\<l 
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Clearly, then Rn = 0(N 1 ' 2 ), and thus it becomes negligible compared to 
the local equilibrium part when N — > oo. 

5.2. The unpinned self-consistent chain. The unpinned self-consistent 
chain with v = was not considered in [11], mainly because the decay of 
correlations is no longer exponential which complicates the analysis. It is 
not completely straightforward to see that the decay is sufficiently strong 
for the previous argument to work, and one has to consider the right subset 
of observables to find strong enough decay; even at equilibrium (qjqj') = 
T((— A) -1 )jj/ which implies that Ylj j'iljlj') = ^(A^ 3 ). However the energy 
only depends on the r x variables, and these turn out to have better decay 
properties. We show in this section that (|5.ip holds also in the self consistent 
profile for v = 0: its first term is 0(1) and the remainder Rn has a bound 
0{N~ 1 / A In 2 N) and is thus dominated by the first, local equilibrium, term. 

Even with v = 0, the formulae for the steady state covariance given 
in Section 2 of [11] hold, and each component is a continuous function of 
v > 0. Thus we can make a shortcut by considering the v — > limits of the 
earlier derived expressions for the steady state expectations. Explicitly, we 
get from Sections 2 and 4 of [H] that for any i,j G {1,2, .. . ,N}, v > 0, 
and any temperature distribution {T n }, there is a unique steady state with 
covariance matrix satisfying 

N N 

(5.22) ( m ) = Y,B[%j)T n , (q iPj ) =^B^(i,j)T n 

n=l n=l 
N 

n=l 

with the explicit definitions of B^"} given in Appendix [Dj We recall the 
definition of r x = q x+ \ — q x given in Section [37171 as well as the representation 

2 2 f 2 

8 X = ^f+-f-\-^p- of the local energies in terms of these variables. Therefore, 
we have 

£j = \v) + Vi^r , with ig = l[x = y]±(l[x = j] + l[x = j - 1]) , 

and thus in order to study the local energy correlations it suffices to con- 
sider the correlations of variables r and p, following the computations in 
the previous subsection. We also point out, that since these variables are 
independent in the equilibrium Gibbs measure, we have (Sf^j)^) = 
unless \f — j\ < 1. 

The main new ingredient allowing us to extend the result of [11] to the 
unpinned case are the bounds proven in Appendix [Dl There we show that 
the exponential decay of the pinned correlations will be replaced by a pow- 
erlaw decay, whose strength depends on the observable. Explicitly, we prove 
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there that in addition to (I5.22D also 



A 1 



N 



(5.23) (r x r y ) = Y J Bl n '(x,y)T n , (r xPj ) = £ B$(x,j)T n , 

n=l n=l 

and the relevant S-matrices decay in m := 1 + \x — y\ + \x — n\ + \y — n\ so 
that there is a constant cq independent of L, x, y, n such that for v = 



(5.24) \B$(z,y) 



< c m 2 (1 + lnm) , B^l{x, y) < com s (l + lnm) , 



< 



com 



Instead of going through all the details here, let us only list the changes 
necessary in the argument used in [11]. As mentioned before, for any tem- 
perature profile there is a unique Gaussian steady state, and the explicit 
formulae for its covariance matrix are still valid for v = 0. In addition, 
Section 3 of [11] also holds verbatim, since the key bound in Lemma 3.2 
does not use the pinning. This shows that also in the unpinned case there is 
a unique self-consistent temperature profile which is bounded between the 
fixed boundary temperatures. The main changes needed will be to Section 
4, since the above bounds produce an additional logarithmic term in the lo- 
cal equilibrium approximations. By the bounds in (|5,24p we have for v = 
and a = 3, 5, 6 that there is a constant d such that 



(5.25) 



N 

E 

n=l 



B^(x, y )\<c'(l + \x-y\) 



3 

4 < 



1, 

(1 + \x 
(1 + Is 



y\ 
y\ 



(the power is not optimal but will be sufficient here) and 

l 1 

(5.26) 



n 



n=l 



B^(x,y)\<c' In 2 N{(l + \x- y \)-\ 

(l + \x-y\)~ 2 , 



if a = 5 , 
if a = 6 , 
if a = 3 , 

if a = 5 , 
if a = 6 , 
if a = 3 . 



Thus, in particular, for X X ,Y X S {r x ,Px} the equilibrium and the self- 
consistent expectations satisfy 



(X x Y y )^ <CT(l + \x-y\) 
\{X x Y y )\ <Cmax(^,T r )(l + |x 



V 



-b 



where 6 = 0, 1,2, depending on the choice of X and Y. This proves that 
{£ji\£j)^' T > is uniformly bounded in N, and since it is zero unless j\ < 1 
we can conclude that (|5.6|) holds in the selfconsistent profile also for v = 0. 



Finally, denoting := max n \T n+ \ — T n \ we have 
(5.27) 

(X x Y y ) - (X x Y y )^ T ^ < Ce N {l + \x- y\)-\l + In 2 iV + \x - j\ 
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We can now use these bounds to complete the proof of Fourier's law and 
derive an estimate of £n for the selfconsistent profile, following section 4.3 
of [UJ. There it was shown that the total current only depends on the 
expectation (qiqi) — {qnQn)- Since this is equal to (roro) — (r^vrjv), we can 
conclude from the above estimates that the discussion on page 795 of [TT] 
holds verbatim, and therefore also for v = we have £n = 0(N~ l l 2 ) in the 
selfconsistent profile. (There will be an additional In 2 N term multiplying 
the right hand side of (4.15) in [TTJ, but this will not change the conclusions.) 

We can now follow the same steps as in Section 15.11 and conclude that the 
bound in (|5.12p for the remainder holds, after we change A to A and q to r on 
its right hand side. Then the above estimates can be applied there, proving 
that the term is bounded by TV -1 / 4 In 2 N. The worst term is the fourth one 
which only has decay \x' — x| -3 / 4 arising from the NESS expectation (r x /r x ), 
and thus leads to a bound e^N 1 ^ In 2 N. This concludes the proof of the 
claims made in the beginning of this section. 



5.3. Higher dimensions. The easiest way to extend the selfconsistent sys- 
tem to a square lattice in higher dimensions is to label the particles with 
x G {1,2,..., N} d =: 1^, and then use fixed boundary conditions in the 
first direction, as before, and periodic boundary conditions for the harmonic 
dynamics in the other directions. We have shown in [11] that this system has 
a unique selfconsistent profile if all heat baths attached to x with x\ = \ are 
fixed to temperature Tj? and those with x\ = N are fixed to temperature T r . 
Moreover, this profile only depends on xi, and the different Fourier- modes 
in the orthogonal directions decouple. 

Explicitly, if we define for k G 1^ _1 , j G ij\r, 

(5.28) %(k):=N-^ £ q{j , y) e- i2 ^ N , 

and pj (k) similarly, then in the selfconsistent steady state 

(5.29) (q^k'ym)) = <^) (1;k) l(k' = k) , 

where the second mean is taken in the stationary state of the one- dimensional 
selfconsistent chain with pinning £ 2 (k) := v 2 + 4 Yli=i sin 2 (nki/N) > v 2 . 
Analogous equations are valid for (<fj/(k')*<7j(k)) and {pj/(k')*pj(k)). 

Using the above partially Fourier-transformed variables, we can write the 
total energy as 

(5.30) 

iv+i , 
M=EE [|p i (k)| 2 + P 2 (k)|^(k)| 2 + i £ |^(k)-?,(k)| 2 . 
*=° keil- 1 \i'-i\=i 
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Therefore, 

N+l 

(5-31) {H) s , s = E(^( k )) (1;k) ' 

where 

(5.32) ^(k) : =^ +l >(k) 2 |+i £ (^-%) 2 - 

i': |i'-i|=l 

In addition, using the independence of the modes and the property (^(k)* 
<&(— k), we also find 



AT+1 

^ \'"J '-/S,S — jy-rf-l ^ ^ 

keI d-i i>,i= 



(5-33) i 3 ( W; «> M = -^ r £ 1 ^(k);^))^) 



Therefore, the energy fluctuations are given by a convex combination 
of the previous one-dimensional fluctuations, only with a varying pinning 
parameter, z^(k) 2 > u 2 . It follows that the statements made in the beginning 
of this section generalize to higher dimensions, i.e., we can conclude that 

(5-34) _L(^) SiS = _L £ (£ x ,;£ x )^ T i c ) + R N , 

x',xe/^ 

where the correction term to local equilibrium, Rpj, is of order l/y/~N in the 
pinned case and iV -1 / 4 In 2 iV in the unpinned case. The second bound could 
likely be improved by using (|5.33j) , but since neither of the one-dimensional 
bounds is optimal, let us skip such computations. 
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Appendix A. The NESS of the velocity-flip model 

Let us first consider the harmonic chain in contact with only two heat 
baths. Its generator is A + £>i,t,> + Bn,t t - It is easy to show that if we start 
the dynamics from an initial centered Gaussian state /2o then the law of the 
process at time t is given by a centered Gaussian state fit. We shall denote 
by C(t) the 2N x 2N covariance matrix of fit- It can be written in the form 



C(t) 



U(t) Z(t) 
Z*(t) V(t) 



HARMONIC CHAINS WITH BULK NOISES 



23 



where U,V, Z are N x N matrices defined by 

Ui,j(t) = J qiqjdfk, Vij(t) = J PiPjdfrt, Zij(t) = J qiPjdJk- 

Let T = (Ti + T r )/2 and 2r\ = T _1 (Ti — T r ). We introduce the matrices A 
and D defined by 

-l\ (0 

$ R ) ' V 2T{R + r]S) 

where the A x N matrices R, S and $ are given by Rj = + 5i,jv)j 

'S'ij = ~ k,N) and = (u 2 + 2)5ij - 5{ + ij - 

Then C(i) is the solution of the differential equation (see |23| ) 

(A.l) ^C(t) = D - AC - CA*. 

dt 

As t goes to +oo, C(t) converges to a positive definite symmetric matrix 
Coo satisfying D = AC^ + C^A* . We note that C^ = TC CCL (1) + rjC where 
C is independent of Tg,T r . 

We now define a Markov process (C(t))t>o with state space S, the set 
composed of the 2N x 2A^ symmetric non-negative matrices C. This Markov 
process is a kind of dual process of the velocity flip process. For any x € 
{1, . . . , A^} we shall denote by T x the diagonal matrix of S such that (r x )j j = 
1 — 2<5j.7v+a;- Consider A independent Poisson processes Af x , x = 1, . . . , N, 
with rate 7/2 and denote by t x (1) < . . . ,r x (k) < . . . the successive times 
at which M x jumps. Let t x < t X2 < . . . be the sequence composed of all 
these times ordered in increasing order. The notation used is such that T Xj 
is a jump of the Poisson process M Xj , i-e., T Xj = T Xj (kj) for the kj-th jump 
of the site Xj. 

The time evolution of (C(t))t>o is given by the following rules: For t Xj < 
t < Tx^+t the evolution of C(t) is deterministic and prescribed by (jA.lj) . 
At time t = T x f +1 , the value of C(t) is given by r :rj+1 C([ri J 4 f 1 1 ]")r :!;j . +1 . It is 

easy to check that the generator C of (C(t))t>o is given, for any continuously 
differ entiable function F : X — > R, by 

2N 

(CF)(C) =J2(D-AC- CA% y d Cx , y F (C) 
(A.2) x ' y=1 

y ' N 

+ |E [ F (r x cr x ) - f{c)\ . 

x=l 

The law of (C(t))t>o when starting from Co is denoted by pt(Co, dC) and 
the centered Gaussian measure on fi^r with covariance matrix C is denoted 
by G c . 
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Lemma 1. Consider the velocity flip model and assume that the initial state 
is jiQ = Gc ■ Then the law \xt of the process at time t is given 

IH= I G cP t(C ,dC). 
Jt, 

Proof. The dynamics generated by C can be described as follows. Let 
(Nx)xe{i,...,N} be the sequence of independent Poisson processes on (0, +oo) 
introduced above. During the time interval [ri j , rlf +1 ) the process (w(t))t>o 
follows the evolution prescribed by A + Bi t T e + &l,T r and, at time rlf +1 , the 
new configuration is obtained by flipping the momentum of particle at site 
Xj+i. Then the system starts again following the dynamics generated by 
A + Bi t T e + Bi t T r until the time of the next flip, and so on. 

Thus, conditionally to the realization of the (t^-)j>i, starting from Gc , 
the law of oj{t) at time t is Gaussian with a covariance matrix C(t) which is 
obtained by the following scheme: in the time interval [7xj>Txj+i)> C(t) sat- 
isfied (lA.lj) ; at time t x ^ +1 , the covariance is given by T Xj+1 C([T^ 1 ]~)r Xj+1 . 
This follows from the fact that if oj has law Gc then u x has law Gc with 

c = r^cr^. 

Hence, the conditional law of the covariance matrix is given by the law 
of the process (C(i))t>o conditionally to the realization of the Poisson pro- 
cesses. The result follows. □ 

We can now prove the following 

Proposition 3. There exists a probability measure Ps } cov> which is invariant 
for (C(t))t>o and whose support is included in the set £ composed of definite 
positive symmetric matrices, such that is given by 

(■)s,v = J Gc Ps,cov(dC) 

Proof. We assume v > 0. The case v = can be treated similarly by 
considering the new variables rj = qj+i — qj. The existence of an invari- 
ant probability measure for (C(t))t>o is based on the existence of a good 
Lyapunov function. 

For any a > we consider the function W a : Q n — > (0, +oo) defined by 

W a {uj) = exp{a%(u;)} 

It can be proved (see e.g. [22] or [5]) that, if a > is sufficiently small, then 
there exists a constant Kq > such that for any initial condition ojq € O n , 

supE^ [W«(w(t))] < K (W a (uj ) + 1) 

where (ui(t))t>o is the process generated by A + B\ t T t + £>7V,T r - Since a flip 
of p x does not modify the value of the energy H{u), this inequality remains 
in force with oj(t) replaced by the velocity flip system uj(t): 

(A.3) supE^o [W Q (w(t))] < K (W„(o> ) + 1) • 
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By integrating ooq with respect to Gc , we get 

sup / Pt (C , dC) G c [W a ] < K [G Co [W a ] + 1] . 
t>o J 

We choose a matrix Co such that the RHS of the previous inequality is 
finite and we remark that W a > aH so that 

sup / p t (C ,dC)GclU}<K' 
t>o J 

for a suitable positive constant K' . Let us write the matrix C as 

C=( u z 
\ Z* V 

and observe that since C 6 E, maxjj \Cij\ < maxjjc^j, V^j}, Ui^jV^i > 0. 
We have G C (U) = tr(J7$) + tr(V) = tr($ 1 /2 [ /$i/2) + tr (y^ ^hus, there 
exists a constant Kq such that 



sup / p t (C ,dC) 
t>o J 



max \Cij\ 
id 



<K>>. 



This shows that the sequence of probability measures {pt(Co, dC) ; t > 0} 
is tight. Consequently, any limiting point Ps,cov 

of the (tight) sequence 



1 



p t (C , dC) := - / p y (C , dC)dy;t>0 



n 



is an invariant probability measure of the process (C(i))t>o- Moreover, it is 
easy to show that J Ps,cov 

(dC) Gc is equal to (-) s ,v This follows from the 
fact that if / : Q n — > M. is a bounded continuous function then F:CeS-> 
Gc(f) is a bounded continuous function on S, so that, if t' is a sequence of 
times such that p t > converges to /c s ,cov, we have 

(/) s , v = lim ~ f dyp y (f) = lim ^ f dy f p y (C ,dC) G c (f) 

t'^co V Jq t'^oo V Jo J 

= J^PAF) = p s , cov (F) = jT f(u) yj p s , C o V (dC) Gc(du)^j . 

It has been proved in [5] that (-) s ,v has a density with respect to the Lebesgue 
measure. This implies that the support of /0 SjCOV is contained in the set S. □ 



Appendix B. Stochastic calculus 

Here we give a proof of results we use to get in one of the steps involved 
in obtaining hydrodynamical equations below equation (|3.13p and the iden- 
tities for the martingales following equation (14. 4j) We begin by stating the 
following standard lemma whose a proof can be found for example in |19j 
or [21]. 
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Lemma 2. If A(u), B(oj) are two local functions of the configuration u 
(r, p), then 



A(u(tN 2 )) - A(u(0)) = N 2 £(£A)(uj(sN 2 ))ds + M t A 

fo< 



B{oj(tN 2 )) - B(co{0)) = N 2 /'(£B)(w(siV 2 ))ds + M? 



where M^,M t B are two centered martingales with quadratic variations given 
by 



(M t A M t B ) = / ds [Mu z (sN 2 ) - A(u(sN 2 ))] 

x [B(oj z (sN 2 ) - B(oj(sN 2 ))} . 

The fluctuation-dissipation equations give 

£{£x) = A4>x+i + C(h x - h x+ x), 
£(r x ) = Ar x + -/~ 1 C(p x - p x+1 ) 

We apply Lemma [2] with 

A(u) = -^^F(x/iV)(r :c + 7 - 1 Vp :[ .), 
(B.2) i x 

B(u) = —J2G(x/N)(£ x + Vh x ) . 



N- 

Jl 

Let F, G : [0, 1] — > K be two test functions vanishing on the boundary. Then 
(B.3) N-^^Fix/N)^ + 1 - 1 Vp x )(tN 2 ) 

X 

- N-V 2 J2 F (x/N)(r x + "f-'VpxM 

X 

« 7" 1 f N- 1 / 2 V F"(x/N)r x (sN 2 )da + A<f (i) 
Jo 

and 

(B.4) N-y 2 Y,G^/N){S x + Vh x ){tN 2 ) 

X 

- N-^Y^Gix/N^+VKm 

X 

« 7" 1 T iV- 1/2 V G"{x/N)4> x {sN 2 )ds + M^i) , 
Jo 
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where M.f are two martingales whose quadratic variations are given 

by 



(Mf (i)Mf (*)) = ^ / ds [^( S iV 2 ) - ^( S iV 2 ))] 

2 ^ 

(M 2 w (t)Mi v (t)) = VE rfs [B(^ z ( 5 iV 2 ) - ^^(siV 2 ))] 

, Jo 



2 

r2 



(M?{t)M»(t)) = / d S [^(siV 2 ) - ^( S iV 2 ))] 

x [^(^(siV 2 ) - ^(siV 2 ))] 
In fact, the two terms 

N- 1 ' 2 Y, F(x/N)(V Px )(tN 2 ) - N- 1 ' 2 Y F(x/N)(V Px )(0) 

X X 

and 

N ~ 1/2 G(x/N)(Vh x )(tN 2 ) - N" 1 ' 2 Y G{x/N)(Vh x )(0) 

X X 

are small in N (just perform a discrete integration by parts and use the 
smoothness of F and G). Using again the smoothness of the functions F 
and G a simple computation shows that 

(B.5) {[M? (ttf) » ^Yf\ F ') 2 (z/N)(p 2 z (sN 2 ))ds, 

([Mm 2 ) « f Q ds{G'f{z/N) 

x (p z (sN 2 ) (r z ( S N 2 )+r z ^( S N 2 )f) , 
(M?(t)M?(t)) « £ MF')(z/N)(G')(z/N) 

x <p 2 (siV 2 ) (r,(siV 2 ) + rw(siV 2 ))) . 

This justifies the limits for the martingales following equation (j4.4p listed in 
Section I4TTTT1 

Appendix C. Proof of equation (14.61) 

Since as t — > oo, u(q,t) — > and e(q,t) — > T(q), (lZ ss ,y ss ) is the station- 
ary solution of 



d t y = (2 7 yw q y-d q ^w 2j 

We prove the statement in (|4.6p only for the energy fluctuations field, the 
proof of the other equation being similar. 
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Let G(q,t) = Gt(q) be the solution of 



(C.l) 



8 t G = (2 7 )- x ^G, 
G(0,t) = G(l,t) = 0, 
[G(q,0) =G (q). 



A 

Then G t (q) = e^^°G (q), and we have 



1 



(C.2) y t (G ) = y (G t ) + — 



V7 Jo J [o,i] 



f(q)(d q Gt- s )(q)W 2 (dq,ds). 



As t goes to infinity, Gt goes to the zero function so that yo(Gt) goes to 0. 
The second term converges to a Gaussian variable with variance equal to 



1 

7 Jo 

-f 

7 Jo 



£(<?„) = -/ da dqf\q)(d q G s ) 2 {q) 
J[0,1] 

ds [ dq{GM{T\q)(d 2 q G s )(q) + 2T(q)T'(q)(d l} G s )(q))} 
J [0,1] 



-2 / / dqf 2 (q)G a (q){d a G a ){q) 

JO J [0,1] 



~(T r -T £ ) ds dqT{q){d q G s )(q)G s {q) 



[0,1] 



dgf 2 (g){GL(g)-G2(q)} 



[0,1] 



~{T r -Ti) ds dqT(q)(d q G 2 s )(q) 
[0,1] 



[0,1] 



dgr 2 (g)G 2 (g) + -(T r - T t f \ ds \ dqG 2 s {q) 



1 



[0,1] 



S A 

since G s = e 2 ^ °G , it follows that 



ds I dqG 2 s {q) 

'[0,1] 



/ d g ( e ^ G )(g) (e^ Ao G )(g) 

'[0,1] 



ds / dq[e-*°G ) (q)G (q) 
'[0,1] 



7 / dgG ( g )((-A )- 1 Go)(g). 

J[0,1] 



This completes the proof of second equality in (|4.6p . 
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Appendix D. Decay of correlations in the unpinned 

SELFCONSISTENT MODEL 

We begin from the expressions given in (I5.22|) . By the results derived in 
Section 4 in [11], we have there for a = 1,2,3 and v > 

(D.l) 

= 9(i -n,j-ri) + g(i + n,j + n) - g(i -n,j + n)- g(i + n,j - n) , 
with g = gN,a,u depending on L = 2(N + 1) and given by 
(D.2) g{ri,n):= ^ f a , u {n' + m' L,n + mL) . 

m' ,m£Z 

The function is obviously L-periodic, so may we assume that its arguments 
have been modified so that \n\, \n'\ < N + 1. Here / denotes a discrete 
Fourier transform of an analytic function whose explicit form depends on 
the "block" a: for all x, y 6 Z, we define 

(D.3) fa,v&,x) := ^ jT dp'£ dpf a Ap',pW {p ' x ' +px) 
with z = 1 — cosp, to = 1 — cosp', and 
(D.4) /i,„(p',p) 



(D.5) / 2 ,„(p',p) := 



(z — U>) 2 + Z + W + V 2 

z — w 
(z — w) 2 + Z + W + V 2 
(z - w) 2 



(d.6) h, v (p',p) -i-r- 

[z — W) z + z + w + v z 

Therefore, for all x, y € {0, 1, . . . , N}, j G {1, 2, . . . , iV}, we also have (recall 
that r x = q x+1 - q x ) 

N N 
(D.7) (r x r y ) = Y,B^(x, y )T n , (r xPj ) = £ bQ(x, j)T n , 

71=1 71=1 

where B^ and 2?g n J are defined as above, using 
(D.8) / 5 ,,(p , ,p):=(e ip ' -l)(e ip -l) 



(z — t«) 2 + Z + W + V 2 

z — w 



(D.9) := ( e V _ i) 

(z — wj z + Z + 10 + i/ z 

(These formulae hold for x, y S {0, iV} since any of the previous functions g 
is 2(iV+l)-periodic and symmetric under x —> —x, separately in both of its 

arguments, which implies that B^(i,j) and B^ihj) vanish if either i or 
j belongs to {0,iV + 1}.) 

We only need to consider < v < L~ l < 1/4, since v — > is to be taken 
before L — > oo. In all of the following, we assume this to hold. Since for 
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all \p\ < 7r we have | sin(p/2)| > |p|/tt, it follows that z > 2iT~ 2 p 2 , and thus 
(z — w) 2 + z + w + v 2 > 2-ir~ 2 (p 2 + (p 1 ) 2 + z/q), vq = on the whole 

integration region. Similarly, as | sinx| < \x\ for all x, \z—w\ < (p 2 + (p') 2 )/2. 
We also have for p = a+i/3, (3 > 0, a G R, that |e ip -l| = |e _/3 — e~ iQ! | < /3+|a| 
and thus also |e _ip — 1| < e^(/3 + |a|). Such bounds will be used frequently 
and without further mention in the following. 

The worst decay estimate will be obtained when / = fs jV . In this case, it 
is likely that for v = the sum in (ID.2P is no longer absolutely convergent, 
which will complicate the analysis. Since then g(n,n') = g(n',n) it will be 
sufflcent to consider the case \n\ < \n'\ < N + 1, which we shall do in the 
following. We begin by inserting the identity 1 = l[|m'|=|?n|] + l[|m / |<|m|] + 
l[|m'|>|m|], which leads to the following decomposition 

oo 

(D.10) g(n ,n) = f(n ', n) + f{n + a'mL, n + omL) 

cr,cr'e{±l} m=l 
oo M 

+ E E E (f(n +a'(M + l)L,n + mL) + 

cr'e{±l} M=0m=-M 

f(n + a'(M + 1)L, n' + mL)) . 



By construction, in the last sum the first argument has always a larger 
magnitude than the second: since \m\ < M, we have \n' + a'(M + 1)L| > 
ML + L — \n'\ > \m\L + L — \n'\ > \n + mL| + L — \n'\ — \n\ > \n + mL|. 

Let us then estimate the first term, f(n',n), which will turn out to be 
dominant. We denote a' = sign(n'), a = sign(n), k' = \n'\, and k = \n\, and 
recall that then k < k' . By a change of variables p — > ap, p' — > a'p' , which 
leave z and w invariant, we thus have 



f(n',n) 



dp'dp cifte+fcV ) (e^V-l)(e^-l) 
^(2^)2 (z-w) 2 + z + w + v 2 



e 

bl,b'l< 



Computing the roots of the polynomial in the denominator yields 
(D.ll) (w — z) 2 + w + z + u 2 = (w — w_)(w — w + ), where 



1 1 AT 

- ± -VI 



w+ = z ±-Vi-8z- Av 2 . 

2 2 



Here y^z denotes the principal branch of the complex logarithm, with values 
such that Re y/z > and \f—\ = i. Then w + is always at least a distance v 2 
away from [0, 2] and w_ always at least a distance | away; more precisely 

(1) If < z < v 2 , then -1 < W- < -| and -6z^ 2 < w+ < -z^ 2 . 

(2) If v 2 < z < | - ^, then -1 < w_ < -| and -Az < w+ < -z. 

(3) If \ - ^ < z < i, then Reu; ± < -\. 

(4) If z > 1, then |Imw±| > \. 
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For z < j£ < ^ — we define z = max(z,i/ 2 ), which implies that then 
— w+ > z. Then if p' = a + i/3 with /3 = \J z/2, we have /3 < 1, and thus 
also cosh/3 < 1 + /3 2 < 1 + z/2. Therefore, then there are pure constants 
C',C that |0 - z) 2 + w + z + v 2 \ > C'(l -cosa + z/2) > C(a 2 + f3 2 ). By 
possibly adjusting constants, this inequality continues to hold for all z and 
a, after we choose z = zq for z > where zq > is a suitably small but 
fixed constant. Then there is a constant cq > such that f3 > co|p| for all 
\p\ < 7T- Therefore, 

J3 M+ff 



|/V,n)|<Ci/ dp|p|e-^ Tdae 



2 



2 



<c 2 r 'dppe-^' r'\^. 

7o Jo 1 + s 

/*7T 

< C 3 / dppe- Copfc '(l + I11P- 1 ) 

<C 4 (l + fc / )- 2 (l + ln(l + fc / ))- 

This proves that the contribution of the first term has an appropriate upper 
bound. Since both of the arguments in the sum in the second term have 
absolute values between (m — and (m + \)L, this bound also proves 
that the second term is bounded by 

C'J,- 2 {\ + InL) < Ci(l + |n'|)- 2 (l + ln(l + |n'|)) , 

which is also of the appropriate form. 

Let us then consider the first term in the final sum in (jD.lOp . We begin 
the analysis by summing over m, which yields the Dirichlet kernel in the 
integrand, amounting to changing inside the defining integral 

e"*-> £ e i(n+mL)P = ^D M (L P ) , where D M {x) == . 

m=— M Sm ^/ 2 ) 



Since Dm(x) is an even function, this yields the following explicit integral 
to be considered 

dp' dp 

1 



/ 

J\v\ 



2 D M (Lp) [cos((n + l)p) - cos (np)] 



x e 



i(ML+L+cr'n')p'. 



(2 — w;) 2 + Z + U> + V 2 

As |Djvf(ip)| < 2M + 1 and |cos((rc + l)p) - cos(np)| < (1 + |n|)p 2 , by 
using the above deformation of the p' integration contour, we obtain that it 
is bounded by 

C 5 M {I + \n\)(l + k')~ 3 (I + Ink'), 

where k' = ML + L + a'n' > L(M + 1/2). As the sum of ln(l + M)/(l + M) 2 
over M = 0, 1, ... is finite and |n| < L/2, this proves the the first term in 
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the final sum in (iDTOj) is bounded by 

C 6 L- 2 (l + lnL) < C 6 (l + |n'|)- 2 (l + ]n(l + |n'|)). 

Since the only property of n, n' used above was that they are bounded by 
L/2, the same bound is valid also for the second term in the final sum. 

The function r i— > r 2 (l + lnr -1 ) is increasing for r £ [0, 1] and thus the 
above bound is decreasing in \n'\. Since max(|n|, \n'\) > (|n| + |n'|)/2, we 
can conclude that there is a constant C5 such that for all n, n' E Z we have 

(D.12) \g(n ,n)\ <c 5 M~ 2 (I + In M), 

where M = 1 + |n mod L\ + \n' mod L\. This is the analogue of the ex- 
ponential bound quoted in (15. 14[) . Therefore, applying this bound and the 
inequality in (I5.15P to (ID.ip proves the first of the bounds stated in (|5.24l) . 
namely that there is C5 such that 

(D.13) \B$(i,j)\ < c' 5 m- 2 (l + lnm) , 

where m = 1 + \n — i\ + |n — j\ + \i — j\. 

For the other two terms (a = 3, 6) the analysis is slightly easier, since 
the bounds will be summable over the Fourier indices. Consider the integral 
defining f a ,u{n,' n), and swap the integration signs as before. If a = 6 and 
\n'\ > \n\, we decompose the integrand first into regions with z > ^ and 
z < jq- For z < jq we decompose further using (jD.lip and 
z — w + — (w — w + ), yielding 

z — w z — w + 1 



(w - z) 2 + w + z + v 2 Vl - 8z - 4z/ 2 w - w + 
1 / z-w + 



w-w- \ y/i-Sz- 4u 2 

Here for the second piece, as well as for the case z > j^, we can then shift 
the integration contour of p' using f3 = /3q which is a small fixed constant. 
These terms have an exponential bound Ce~^°' n The remaining piece we 
estimate using the same /3 = (5{z) as before, which proves that 

|iUn',n)| < C L-W\ + / dp(z + ^)e~^'\ f dae^^l^ 
< C (e-P°\ n '\ + v z lnz^-^-^'l + (1 + |n'|)- 3 (l + ln(l + |n'|))) , 

where we have used the facts that for \p\ < ^ we have \z — u> + | < C{z + v 2 ), 
and that j3 = v and z < v 2 if \p\ <v. If a = 6 and \n\ > \n'\, we need to swap 
the order of integration, but we can still apply the same decompositions as 
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above yielding 



16 

n 



\hAn,n)\<C[e-M n \+ [ dp (z + v 2 )\p\e~ *N f da -^—^ 
< c L-toW + j,3 e -Hn| + + 



Therefore, \h, v {n\n)\ < C[e~P ok ' + v 3 In v~ l e~ uk ' + (1 + A;')^ 3 (l + ln(l + k'))] 
where k' = max(|n|, \n'\), and thus 

\g(n, n)\ < C[M~ 3 (1 + In M) + e"*^ 2 + L~ 2 zy In u' 1 + L~ 3 (l + In L)] 

where M = 1 + \n mod L| + \n' mod L|, as before. Taking here the limit 
v — > and then following the same steps as before proves that a constant 
for the second of the bounds (a = 6) stated in (|5.24j) can be found. 

For the final remaining case a = 3 we can again use symmetry and assume 
from the beginning that \n'\ > \n\. We decompose then for z < jq 

(z - w+) 2 1 



(D.14) hJp',P) 



Vl - 8z - Av 2 w - w A 



( (z — wa ^ 



+ 7= = = 9 + 2(2 - w+) +w + -w. 

and apply the rest of the argument as in the previous case. This yields a 
bound \f 3 ,„(n',n)\ < C[e~^ k ' \v^~ vk ' + (l+/k')" 4 ] where k' = max(H, \n'\), 
and thus 

\g(n',n)\ < C[M~* + e"* M / 2 + L~ 2 v 2 + L -4 ] 

where M = 1 + |n mod L\ + |n' mod L\. As above, this implies the final 
bound in (|5.24p . for a = 3. 
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